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Abstract 

Using exact results, we determine the complex-temperature phase diagrams of the 2D Ising 
model on three regular heteropolygonal lattices, (3 • 6 • 3 • 6) (kagome), (3 • 12 2 ), and (4 • 8 2 ) 
(bathroom tile), where the notation denotes the regular n-sided polygons adjacent to each 
vertex. We also work out the exact complex-temperature singularities of the spontaneous 
magnetisation. A comparison with the properties on the square, triangular, and hexagonal 
lattices is given. In particular, we find the first case where, even for isotropic spin-spin 
exchange couplings, the nontrivial non-analyticities of the free energy of the Ising model lie 
in a two-dimensional, rather than one-dimensional, algebraic variety in the z = e~ 2K plane. 
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1 Introduction 



The Ising model has long served as a prototype of a statistical mechanical system which 
undergoes a phase transition with associated spontaneous symmetry breaking and long range 
order. The Ising model on a dimension d = 2 lattice has the great appeal that (for the spin 
1/2 case with nearest-neighbor interactions) it is exactly solvable; for the square lattice, in 
the absence of an external magnetic field H, the free energy was first calculated by Onsager 
|U, and a closed-form expression for the spontaneous magnetisation was first derived by 
Yang p|. Later, the zero-field model was also solved for the triangular and hexagonal (= 
honeycomb) lattices; for a comparative review of the solutions on these three lattices, see Ref. 
0. However, in addition to the square, triangular, and hexagonal lattices, which involve 
tilings of the plane by one type of regular polygon, there are also other 2D lattices which have 
the property that all vertices are equivalent and all links (bonds) are of equal length, but are 
comprised of tilings by more than one type of regular polygons. We shall denote the regular 
2D lattices which are comprised of tilings by one type of polygon as "homopolygonal" and 
those involving tilings by more than one type of polygon as "heteropolygonal" . Indeed, as 
will be discussed further below, the (zero-field) free energy and spontaneous magnetisation 
for the Ising model have been calculated for three heteropolygonal 2D lattices. The results 
yield valuable insights into how the properties of the solution depend on the tiling of the 
plane and, in particular, how they change when this tiling involves more than one kind of 
polygon. 

In the present paper, we shall determine the complex-temperature phase diagrams and 
singularities of the spontaneous magnetisation for the three particular heteropolygonal lat- 
tices for which the Ising model has been solved. There are several reasons for studying 
the properties of statistical mechanical systems such as spin models with the temperature 
variable generalised to take on complex values. First, one can understand more deeply the 
behaviour of various thermodynamic quantities by seeing how they behave as analytic func- 
tions of complex temperature. Second, one can see how the physical phases of a given model 
generalise to regions in appropriate complex-temperature variables. Third, a knowledge of 
the complex-temperature singularities of quantities which have not been calculated exactly 
helps in the search for exact, closed- form expressions for these quantities. This applies, in 
particular, to the susceptibility of the 2D Ising model. Complex-temperature properties of 
the Ising model (zeros of the partition function and singularities of thermodynamic and re- 
sponse functions) have been explored in a number of papers for homopolygonal 2D lattices 
and for various 3D lattices §, | g, 0, g, |, 0, fnj [12], HI 0- ( Refs - 1 and P art of Refs - 
P, || deal with the d = 3 Ising model.) However, to our knowledge, analogous investigations 
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of complex-temperature properties have not been reported for heteropolygonal 2D lattices. 

To begin, we need to describe the heteropolygonal lattices. We shall use the standard 
mathematical notation for general lattices |15| . An important virtue of this notation is that 
the mathematical symbol explicitly contains sufficient information for one to reconstruct 
the lattice. A regular tiling of the plane is defined as one involving one or more regular 
polygons, i.e. polygons each of whose sides are of equal length. An Archimedean lattice 
is defined as a regular tiling of the plane in which all vertices are equivalent. Thus each 
vertex has the same coordination number, which we shall denote as q. From the definition, 
it follows immediately that all of the links on an Archimedean lattice are of equal length. 
For this reason, it is natural to consider the case of equal spin-spin couplings, and we shall 
do this here. However, it should be noted that the links on an Archimedean lattice are not, 
in general, all equivalent. For the three homopolygonal Archimedean lattices, it is obvious 
that all links are equivalent, but for the heteropolygonal lattices, this equivalence only holds 
if (a) the lattice consists of only two types of polygons and (b) each edge of the first type 
of polygon is also an edge of the second type of polygon. Among the eight heteropolygonal 
Archimedean lattices, this condition is met only for the (3 • 6 • 3 • 6) (kagome) lattice (see 
below). The mathematical definition of, and notation for, an Archimedean lattice A are 
specified by 

A=cn Pi ) (i.i) 

where p denotes a regular p-gon. The meaning of eq. ( |1 . 1| ) is that as one makes a small 
circuit around any vertex, one traverses first the polygon pi, next P2, and so forth, finally 
traversing p q . Because the starting point is irrelevant, the symbol is invariant under cyclic 
permutations: 

(Pi " P2 ■ • • Pq) = (P2 • P3 ■ • 'Pq • Pi) (1.2) 

and so forth for other cyclic permutations. If a p-sided polygon occurs £ times sequentially 
in this product, one signifies this by p . There are 11 such (d = 2) Archimedean lattices 
[ 15| . Of the 11 Archimedean lattices, three are homopolygonal. These are the well-known 
square, triangular, and hexagonal lattices, which are denoted, respectively, (4 4 ), (3 6 ), and 
(6 3 ). Note that the dual of a homopolygonal lattices (p e ) is (£ p ), so that it is self-dual if 
and only if p = £, which occurs only for p = 4. The other eight Archimedean lattices are 
heteropolygonal. The full set of these lattices is listed in Table 1. (One lattice, (3 4 -6), occurs 
in two enantiomorphic (chiral) forms, which are considered together in this counting.) 

The three heteropolygonal Archimedean lattices for which we shall determine the complex- 
temperature phase diagrams and singularities of the magnetisation are those for which the 
Ising model has been solved exactly, namely the (i) (3 • 6 • 3 • 6), (ii) (3 • 12 2 ), and (iii) (4 • 8 2 ) 
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lattice 


name 


Q 


bip. 


sym. 


(3*) 


triangular 


6 


N 


z — > — z 


(4 4 ) 


square 


4 


Y 


z — > — z, z — > l/z 


(6 3 ) 


hexagonal 


3 


Y 


z -> l/z 


(3 • 6 • 3 • 6) 


kagome 


4 


N 


z — z 


(3 • 12 2 ) 




3 


N 


n 


(4 • 8 2 ) 


bathroom tile 


3 


Y 


2 -> 1/Z 


(3 4 • 6) 




5 


N 


n 


(3 3 • 4 2 ) 




5 


N 


n 


(3 2 • 4 • 3 • 4) 




5 


N 


n 


(3 • 4 • 6 • 4) 




4 


N 


z — > — z 


(4 • 6 • 12) 




3 


Y 


z -> l/z 



Table 1: Table of the 11 Archimedean lattices. The first three are the well-known ho- 
mopolygonal ones, and the remaining eight are heteropolygonal. The entries in the column 
denoted bip. indicate whether (Y,N) the lattice is bipartite. The entries in the column 
denoted "sym." indicate whether the complex-temperature phase diagram and zeros of the 
partition function for the (spin 1/2, zero- field) Ising model have the symmetries z — > — z 
and/or z — > l/z, or neither, denoted by "n". 
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lattices. For the reader's convenience, these lattices are shown in Figs. 1-3, respectively. 
The (3 • 6 ■ 3 • 6) lattice is often called the kagome lattice, after a Japanese word for a similar 
type of basket weave. The (4 • 8 2 ) lattice is often called the bathroom tile lattice. In the 
literature, the (3 • 12 2 ) and (4 • 8 2 ) lattices are frequently denoted by the shorthand symbols 
3 — 12 and 4 — 8. As was observed by Utiyama Hl6| , one way to construct these lattices is 
to start from a checkerboard lattice, replace each "black" square by a square with n v addi- 
tional vertical bonds, and then take certain of the spin-spin couplings either to or oo. One 
chooses n v = 1 to construct the (3 • 6 • 3 • 6) and (4 • 8 2 ) lattices, and n v = 3 to construct the 
(3 • 12 2 ) lattice. Using this observation, Utiyama discussed ingredients for an exact solution 
for the (zero-field) free energy for these lattices |fT6| . These lattices can also be obtained by 
(possibly iterated) decoration and/or star-triangle operations from the three homopolygonal 
lattices [|i~7j 

We recall that abstractly, the duality mapping for a d- dimensional lattice maps fc-cells to 
d—k cells. In our present case, with d = 2, this mapping interchanges 0-cells (vertices) and 2- 
cells (polygonal faces) while taking 1-cells (links) to other links. The duals of Archimedean 
lattices are Archimedean if and only if they are homopolygonal; for the heteropolygonal 
Archimedean lattices, the dual lattices are of a different type, called Laves lattices in the 
mathematical literature |15|, [Tj|. A Laves lattice is defined as a regular tiling of the plane 
involving a single (in general, non-regular) polygon. Thus, the links of a Laves lattice are 
not, in general, of equal length, and the vertices (sites) are not, in general, equivalent. There 
is a one-to-one correspondence, namely that of duality, between the 11 Archimedean lattices 
and the 11 Laves lattices. Clearly, the number of different types of 0-cells (vertices) on the 
Laves lattice A* which is the dual of the Archimedean lattice A is equal to the number of 
different 2-cells (polygons) on A. Similarly, the fact that the Laves lattice consists of only one 
type of (generally non-regular) 2-cell (polygon) follows, by duality, from the fact that its dual 
Archimedean lattice consists of only one type of 0-cell (vertex). The standard mathematical 
notation for a Laves lattice A is 

A=[(fU')] (1-3) 

i=i 

the meaning of which is as follows: on a given polygon, as one makes a circuit around 
its periphery, one passes first a vertex with coordination number qi, next a vertex with 
coordination number q 2 , and so forth, finally coming to a vertex with coordination number 
q p , where p denotes the number of (generally non-equal) sides of the polygon, which is, of 
course, equal to the number of vertices of this polygon. The Laves lattice [(ITj=i <lj)] is 
the dual of the Archimedean lattice (]Tj=iPj) for p = q and pj = qj, j = I. ..p. As an 
example, the Laves lattice [4 • 8 2 ] (commonly called the union jack lattice) is shown in Fig. 
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4 and is evidently the dual of the Archimedean lattice (4 • 8 2 ) (bathroom tile). From the 
duality relation with the Archimedean lattices, or directly from the definiton, it follows that 
there are three Laves lattices for which the vertices are, in fact, all equivalent; these are 
[3 6 ], [4 4 ], and [6 3 ]. We shall denote these as homovertitial, and the remaining eight Laves 
lattices with inequivalent vertices as heterovertitial, by analogy with the homopolygonal and 
heteropolygonal Archimedean lattices. 

The solution for the (zero-field) free energy on each of the Archimedean lattices immedi- 
ately gives the solution on the corresponding dual Laves lattice. Thus, our determination of 
the complex-temperature phase diagrams for the three heteropolygonal Archimedean lattices 
considered here also yields a corresponding determination of the complex-temperature phase 
diagrams for their dual Laves lattices. 



2 Model and Notation 

Our notation is standard, so we review it here only briefly. We consider the spin 1/2 Ising 
model on the 2D lattices A as discussed above at a temperature T and external magnetic 
field H (where H = unless otherwise specified) defined by the partition function 

Z =Y j e- m (2.1) 

{a„} 

with the Hamiltonian 

h = -j £ <w-i?5> n (2.2) 

<nri> n 

where a n = ±1 are the Z 2 spin variables on each site n of the lattice A, (3 = (ksT)^ 1 , J is the 
exchange constant, < nn' > denote nearest-neighbour sites, and the units are defined such 
that the magnetic moment which would multiply the H J2 n a n is unity. (Hereafter, we shall 
use the term "Ising model" to denote the spin 1/2 Ising model unless otherwise indicated.) 
We use the usual notation K = (5 J, h = j3H, v = tanhif, z = e~ 2K , u = z 2 = e~ 4K , and 
fi = e~ 2h . Note that v and z are related by the bilinear conformal transformation 

l-v 



1 + v 



(2.3) 



We record the symmetries 

K -> -K {v -> -v , z -> l/z , u — > l/u] (2.4) 
It will also be useful to introduce the common abbreviations 

C = cosh(2ir) (2.5) 
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S = smh(2ir) (2.6) 

The reduced free energy per site is / = — j3F = linijv.,-^00 iV" 1 In Z in the thermodynamic 
limit, where N s denotes the number of sites on the lattice. 

3 Some Basic Properties 

We begin by discussing the phase boundaries of the model as a function of complex tem- 
perature, i.e. the locus of points across which the free energy is non-analytic. The physical 
phases of the model include the phase where the Z 2 symmetry is realised explicitly, viz., the 
paramagnetic (PM) phase, and also, for dimension greater than the lower critical dimension- 
ality (i.e., d > 2, for integer d), a phase where the Zi symmetry is spontaneously broken with 
long-range ferromagnetic (FM) order, i.e. a nonzero spontaneous (uniform) magnetisation, 
M. For the lattices considered here which are bipartite and hence involve no frustration 
for antiferromagnetic (AFM) ordering, there is also a phase with AFM long range order, 
i.e. a nonzero staggered magnetisation, M st . As one can see from Table 1, of the three 
heteropolygonal lattices considered here, only the (4 • 8 2 ) lattice is bipartite. One defines the 
complex-temperature extensions of the physical phases by analytic continuation in K or an 
equivalent variable such as v, z, or u. There are, in general, also complex-temperature phases 
which have no overlap with any physical phase. As in our earlier work, we label these as O 
phases (where O denotes "other"), including subscripts to distinguish them where there are 
several. In cases where one O phase is precisely the complex conjugate of another, we label 
the one with Im(z) > by O (with subscript where necessary) and its complex conjugate 
by O*. (From eq. ( |2.$ ), h follows that in the v plane, the corresonding O and O* phases 
have Im(v) < and Im(v) > 0, respectively.) 

As noted in Ref. ||, there is an infinite periodicity in complex K under the shift K — > 
K + niTT, where n is an integer, and, for lattices with even coordination number q, also the 
shift K — > (2n + l)z7r/2, as a consequence of the fact that the spin-spin interaction a n a n > in 
"H is an integer. In particular, there is an infinite repetition of phases as functions of complex 
K. These repeated phases are reduced to a single set by using the variables v, z or u, owing 
to the symmetry relation K — ► K + nm =^ {y — > v , z —>■ z, u — > u}. It is thus convenient 
to use these variables here. 



As we have discussed in our earlier works [TT|, [T2|, [T3|], the equations for the locus of points 



where the free energy is non-analytic also serve to define the boundaries of the complex- 
temperature phases of the model. Some of the loci of points where the free energy is non- 
analytic do not actually separate any phases, but rather are arcs or line segments protruding 
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into various phases. In passing we note that the free energy is, of course, also trivially non- 
analytic at K — ±00, i.e. v = ±1 or z = 0, 00. This is obvious from eqs. fl2.1|) and (|2.2p . 
Because these points are isolated, they do not separate any complex-temperature phases and 
hence will not be important here. 

Before proceeding, we list some theorems which formalise results in our earlier work. 
Although these are elementary, it will be useful to record them here for our later discussion 
of the complex-temperature phase diagrams. Some of the theorems will be given in greater 
generality than is needed here. 

Theorem 1 

Consider the (spin 1/2, zero-field) Ising model on an arbitrary lattice. The loci of points 
in the complex z and v planes where the free energy of the Ising model is non-analytic 
are symmetric under complex conjugation, i.e. under reflection about the Im(z) = and 
Im(v) = axes, respectively. Hence the same is true of the complex-temperature phase 
boundaries in these variables. The same is also true of the zeros of the partition function on 
finite lattices. 

Proof 

One can prove the first two parts of this theorem either by starting with finite lattices 
or by considering only the properties of the thermodynamic limit. If one uses finite lattices, 
this automatically also proves the last part of the theorem, concerning zeros of the partition 
function. The basic property that one uses is that the partition function on a finite lattice 
is a generalised polynomial. (Here, a "generalised polynomial" in the variable ( is defined as 
function consisting of a finite sum of integral powers of (, where negative as well as positive 
powers are allowed. Specifically, the partition function has the form 

Nt/2 

Z= £ c pZ P (3.1) 

p=-N e /2 

where Ne denotes the number of links on the lattice and the c p are real. (In fact, c p are not 
just real, but integral, and also satisfy the symmetry c p = cn ( - p , but these properties will 
not be needed here.) The zeros of Z are the zeros of the polynomial 

P( z ) = z N ^z = c P z p (3.2) 

p=0 
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Since the polynomial P{z) has real coefficients, it satisfies 

P(z)* = P(z*) (3.3) 

Hence, if Zj is a zero of P(z), then so is z*. This proves the third part of the theorem, 
on the zeros of the partition function for finite lattices. Taking the thermodynamic limit, 
one obtains the result that the non-analyticities of the free energy and hence the complex- 
temperature phase boundaries, are invariant under z — > z*. From the mapping fl2.3|) relating 
z and v, the same result follows for the non-analyticities and phase boundaries in the v plane. 
□ 

In passing, we note that if one uses periodic boundary conditions, then the number of 
links is simply given by Ni = (q/2)N s . Moreover, for the part of the theorem which deals 
with the free energy (in the thermodynamic limit), an alternate way of seeing the result is 
to observe that / has the form 

/ = ln2 + 141 + «(*,«,)] (3.4) 

where A and B are rational functions of z, or equivalently, v, with real coefficients. Aside 
from isolated points at K — ±00 (i.e., v = ±1, z = 0, 00) and, for odd q, the point z = — 1, 
the above locus of points is the set where A + 2 ) — 0. But since this amounts to a 
polynomial equation in z or v with real coefficients, it follows that the roots of the equation 
are either real or come in complex conjugate pairs. 

This theorem can be generalised in two ways. Although we shall not need either of these 
generalisations here, we state them for completeness. First, the theorem actually applies not 
just for the case of zero external field, but also for the case of nonzero (real) magnetic field. 
One can also show that it is true for pure imaginary h. Secondly, one can straightforwardly 
generalise Theorem 1 to the case of the spin s Ising model. For this purpose, we define the 
normalisation of the spin values to be a n G {—2s, — 2s+2, 2s— 2, 2s} in eqs. ( |2.1| ) and (|2.2| ). 
Then the generalisation of Theorem 1 reads as given there, without the restriction to spin 
1/2. For integral s, a commonly used alternate normalisation is cr n G {— s,— s + l,...,s — l,s}; 
with this normalisation, the generalisation of Theorem 1 would use the modified definition 
z = e- K .) 

Theorem 2 

For the (spin 1/2, zero-field) Ising model, if and only if the lattice has even coordination 
number q = 2r, the complex-temperature phase diagram is invariant under the transforma- 
tion z — > —z. The same is true of the zeros of the partition function for finite lattices. 
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Proof 



This follows by explicit calculation of the partition function for finite lattices. For odd q, 
Z is a generalised polynomial in z, while for even q, it is a generalised polynomial in u = z 2 . 
Using the definition of the free energy / = timjv^oo N" 1 In Z, it follows that the locus of 
points where the free energy is non-analytic has the same symmetry. □ 

Thus, for lattices with even q, a more compact way to display the complex-temperature 
phase diagram is in the complex u plane. 

Theorem 3 

For the (spin 1/2, zero-field) Ising model, if and only if the lattice is bipartite, the 
complex-temperature phase diagram is invariant under the transformation z — > 1/z {u — > 1/u 
if also q is even). The same is true of the zeros of the partition function for finite lattices. 

Proof 

The statement that the lattice is bipartite means that it can be decomposed into two 
sublattices, which we may denote as even (e) and odd (o), such that each site on the even 
sublattice has as its nearest neighbors only sites on the odd sublattice, and vice versa. Since 
the spin-spin interaction in the Hamiltonian only involves nearest-neighbor pairs, it follows 
that under the mapping 

<7 e — > <7 e 
a a — > — a a 

J -> - J (3.5) 

the partition function Z and hence also the free energy / are invariant, where a e and a Q 
denote spins on the even and odd sublattices. Hence 

Z(K) = Z(-K), f(K) = f(-K) (3.6) 

The theorem follows, since K — > — K is equivalent to z — > 1/z. □ 

Theorem 4 

In the z plane, the zeros (and any divergences) of the spontaneous magnetisation M 
occur at either real values or at complex conjugate pairs of values. 
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Proof 



Theorem 4 is a corollary of the generalisation of Theorem 1 to the case of nonzero (real) 
external field h. Since M = lim^^o dM/dh, this generalisation of Theorem 1 implies that 

M(z)* = M(z*) (3.7) 

Hence, in particular, the set of zeros of M as a function of z is invariant under z — > z*. 
Similarly, the set of values of z where M diverges (if this set is non-null) is invariant under 
z^z*. □ 

The following result is a corollary of this theorem: 
Theorem 5 

On a bipartite lattice, in the w = 1/z plane, the zeros (and any divergences) of the 
spontaneous staggered magnetisation M st occur at either real values or at complex conjugate 
pairs of values. 

Proof 

This theorem follows from the fact that the symmetry ( |3.5| ) maps the model with the 
ferromagnetic sign of the spin-spin coupling, J > 0, to the model with the antiferromagnetic 
sign, J < 0. This allows one to obtain the staggered magnetisation M st immediately from the 
uniform magnetisation; in the physical AFM phase and, by analytic continuation, thoughout 
the full complex-temperature extension of the AFM phase, M st is given as 

Mjw) = M(z -> w) (3.8) 

where w = 1/ z = e 2K = e" 2 '^"'. (M st vanishes identically elsewhere.) □ 

Theorem 6 

For lattices with odd coordination number q, the zero-field partition function Z of the 
(spin 1/2) Ising model vanishes at z = — 1, and the free energy contains a negatively divergent 
singularity at this point. 

Proof 
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From the definition of the partition function and of K, we have 



Z = exp(K ^ cr n a n >) (3.9) 

{<r n } <nn'> 

Now, in general, for each link < nn' >, since <J n a n > = ±1, it follows that 

e K* n a n , = coshK + 0nGnl sinhK (3.10) 

The K values corresponding to z = lim e ^ ( — 1 =t e ) (where e denotes a small real number) are 
K = —(1/2) lim e _^o ln(— l±e) = — (l/2)(±i7r+2m7r), where we follow the usual convention of 
taking the branch cut for log z to lie along the negative real z axis, and n indexes the Riemann 
sheet of the logarithm. To minimize unimportant minus signs, we consider z = lim e ^o( — 1 — e) 
and take the principal Riemann sheet of the log, n = 0; then K = m/2. Substituting this 
into eq. ( |3.10p , we have 

e i**/2)w = i an(Tnl (3.11) 

so that 

Z = E ( II «w) (3.12) 

{a n } ^<nri> ' 

Since the lattice has coordination number q, in this product over all links, the a n for each 



site n appears q times, and hence eq. (|3.12| ) can be rewritten as a product over all sites: 



^EIK) ( 3 - 13 ) 



Evaluating the summation over each spin, we have 

:N, 



+iy + (-M Ns (3.i4) 



Evidently, for odd q, the summation over o n = ± on each site yields 0, and hence Z = 0, so 
that the free energy is negatively divergent at this point. □ 

Among the lattices considered here, the hexagonal, (3 • 12 2 ), and (4 ■ 8 2 ) ones have odd 
q (in each case, q = 3), so that Theorem 6 implies that the respective free energy / for each 
has a divergent singularity at z — — 1. 

From eq. ( |3.14| ) in the proof of Theorem 6, it is immediately evident that for the (spin 
1/2) Ising model on a lattice of even q = 2r, the value of the partition function on a finite 
lattice with periodic boundary conditions (and hence the relation = (q/2)N 8 ) is 

Z(z = -l;q = 2r)=iT N '2 N ' (3.15) 
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while the (reduced) free energy satisfies 

/(* = -l;g = 2r) = rIni + ln2 (3.16) 

(independent of boundary conditions, since these do not affect the thermodynamic limit). 
In eq. (|3.16|) , \ni = in/ 2 + 2irin, where n labels the Riemann sheet of the log. 

The points across which the free energy is non-analytic are related to the zeros of the 
partition function. A fundamental question concerns whether these points (apart from the 
trivial ones at K = ±00 and, for odd q, at z = —1 as proved in Theorem 6) lie on curves 
(including line segments) in the z (or equivalently, v) plane, or whether they lie in areas. For 
homopolygonal 2D lattices with isotropic spin-spin exchange couplings, these points do lie 
on curves. (There is also numerical evidence for this in the case of 3D lattices ||.) It is also 
well known that for 2D homopolygonal lattices with anisotropic couplings, these points lie, 



in general, in areas rather than on curves [19]. It is worthwhile to investigate this question 



for heteropolygonal lattices, and we shall do so, as part of our general determination of 
the complex-temperature phase diagrams. Indeed, one of our interesting results will be the 
finding that even for isotropic spin-spin couplings, the zeros do not always lie on curves in 
the case of heteropolygonal lattices. We shall give a simple explanation of this finding below. 
We now proceed with our analyses. 



4 (3 • 6 • 3 • 6) Lattice 



The (zero-field) free energy of the Ising model on the (3 • 6 • 3 ■ 6) (kagome) lattice with equal 
spin-spin couplings was first calculated explicitly by Kano and Naya p0|. The result is 



/**. = In 2 + ~ £ £ d -^0 ln{~ [(C 3 + S 5 )" + 3C 2 - 2CS\C + S)P(9 1 , 9 2 )] } (4.1) 
where 

P(0i, 6 2 ) = cos(0i) + cos(0 2 ) + cos(0! + 2 ) (4.2) 

As was the case with the homopolygonal lattices, the connected locus of points across which 
the free energy is non-analytic is the set of points for which the logarithm in the integral 
of eq. ([4.1] ) vanishes. These are also the zeros of the partition function. In terms of the 
variable u, this vanishing condition is the equation 

(21m 4 + 24u 3 + 18u 2 + 1) - 4w(l + it)(l - ufx = (4.3) 
where x represents P(9±, 62) and takes values in the range 

- - < x < 3 (4.4) 
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For x = 3, there is a double real root at u c ^ ag _ } where 

u cMg . = -1 + ~ 0.15470054... (4.5) 

and another double real root at — l/(3it Cj fc 0ff . = — 1 — 2/y / 3- As indicated in the notation, 
u c ,kag. is the physical critical point separating the FM and PM phases. As x decreases from 
3 in the range — 1 < x < 3, each of these splits into two pairs, which trace out the curves in 
Fig. 5(a). For x = — 1, these rejoin in conjugate double roots at the multiple or intersection 
points Uk and u* k , where 

u k = -(-l+2i) = h~ 1/2 e ldk (4.6) 
5 

with 

9 k = n- arctan(2) ~ 116.57° (4.7) 



As in our earlier works [11], [T2|, [13|], we denote these as multiple points, following the technical 
terminology of algebraic geometry, according to which a multiple point of an algebraic curve 
is a point where two or more branches (arcs) of the curve cross . We shall use the words 
"intersection point" and "multiple point" synonymously here. The corresponding numerical 
values of z k = ±5~ 1 / 4 e i6 ' fc / 2 are 

z k = ±(0.3515776 + 0.5688645i) (4.8) 

Finally, as x decreases from —1 toward —3/2, these split again, with one root moving away 
from u k upward along an arc of the circle defined by 

\u-l\ = l (4.9) 
toward the endpoint of this arc on the postive imaginary axis, given by 

Ue = ^ (4.10) 



the conjugate root moving down from u* k on the circle (|4.9| ) toward u* e , and the other two 
roots moving along the circle (|4.9| ) toward the real axis. Finally, at x — —3/2, two roots 
occur at the endpoints u e , w*, and there is a double root at the point u = u e = —1/3. 

In terms of the commonly used variable z, the complex phase diagram is shown in Fig. 
5(b). From the fact that the equation ( f4.3|) for the locus of points where / is non-analytic is 
symmetric under z — > —z, it follows that the phase diagram in the z plane has this symmetry 
also. The physical critical point is z c = (—1 + 2/y / 3) 1 ^ 2 — 0.393320.... The boundary of the 
complex-temperature FM phase crosses the real z axis at the points ±z c . Corresponding to 
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the intersection point u k and its conjugate there are four intersection points z k = 5 1 / 4: e l9k ' 2 , 
— Zk, and rkz^. The outermost points on the boundaries of the O phases are given by 
2o = i|(-l-2/v / 3)| 1/2 and z* . 

Finally, by a conformal mapping or directly from eq. ( f4.1|) , one can also determine the 
corresponding loci of points in the v plane. The equation for this locus of points is 

1 - Av + 10v 2 - 16v 3 + 22v A - 16v 5 + 10v 6 - Av 7 + v 8 - 2v 2 (l - vf{l + v 2 )x = (4.11) 

Up to an overall factor of v ~ 8 , this equation is invariant under the transformation v — > 1/v. 
This together with the reality of the coefficients implies that the locus of solutions is invariant 
under (i) v — > v* and (ii) v — > 1/v. This locus is shown in Fig. 5(c). The complex- 
temperature phases are marked on this figure and consist of the PM and FM phases, together 
with two O phases which are related to each other by complex conjugation. From the value 
of u c or z c and the relation (|2.3| ), it follows that the critical value of v is 

v cMg . = ^(1 + 3 1/2 ) [l - (2^3 - 3) 1/2 ] ~ 0.43542054... (4.12) 

One sees that the kidney-shaped curve enclosing the FM and two O phases crosses the real-v 
axis at the points v c and l/v c ~ 2.29663.... The four intersection points of the kidney- 
shaped curve and the arcs are given by v & = (1 — z k )/{l + z k ), its complex conjugate, t>|, 
and their two reciprocals, 1/vj. and l/v%- Similarly, the four endpoints of the arcs are given 
by v e — (1 — z e )/(l + z e ), v*, l/v e , and 1/v*. We find that the inner endpoints of the arcs 
are the same distance from the origin as the physical critical point, i.e. \v e \ = v c . 

The phase structure consists of PM and FM phases, but, as is well known, no AFM phase 
because of the frustration associated with AFM ordering. These characteristics are the same 
as those for the (isotropic) Ising model on the triangular lattice. In addition, Figs. 5(a) and 
5(b) show two complex-temperature phases denoted O and O* which have no overlap with 
any physical phase. In the u variable (see Fig. 5(c)), these are mapped onto a single O phase. 
In our previous work, we found O phases for the square and triangular (but not hexagonal) 
lattices, in the v and z plane. However, for both the square and triangular lattices, each 
point z in the O phase was just the negative of a corresponding point in the complex PM 
phase, and hence under the mapping from the z to u plane, these were mapped to the same 
point, which could be considered just the complex PM phase. We find for the kagome lattice 
a new property, viz., a distinct O phase which persists even in the u plane. 



The spontaneous magnetisation cannot be written in the same form M = (1 — (k <j 



A) 



\2M/ 



as for the regular unipolygonal lattices A = sq, tri, he. Rather, it has the following form in 
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the FM phase (and vanishes identically elsewhere) 



2 / \ 1 



where 



i.e., 



M kag . = ^ ' + v ^ } ( 1 - (k <Mg f ) (4.13) 



2 7 /V/ 2 (1 + M ) 3 / 2 (1 + 3m 2 ) 1 /2 



(1-6 M -3m 2 ) 1 /8(i + 2m + 5m 2 ) 3 / 8 (1 + 3m) 1 / 4 1 
Mka9 - ~ (l- U )V4(i + u) ( 4 - 15 ) 

These formulas may be analytically continued throught the complex-temperature extension 
of the FM phase. Note the factorizations 1 — 6u — 3u 2 = (1 — u/u c ^ kag ){\ + 3u C:kag u) and 
1 + 2u + 5m 2 = (1 - u/u k )(l - u/u* k ). 

From eq. ( |4.15|) and our determination of the complex-temperature FM phase in which 
the analytic continuation of this formula holds, it follows that besides the well-known fact 
that M kag vanishes continuously at the physical critical point u c , with exponent ft = 1/8, it 
also vanishes continuously at the complex-temperature points ug = —1/3 with exponent 

/W = \ (4-16) 
and at the points u k and u* k with exponent 

/W = I (4-17) 

Elsewhere along the boundary of the complex-temperature FM phase, M vanishes dis- 
continuously. One may observe that the expression for M has an apparent zero at u — 
—l/{3ukag.,c) = — 1 — 2/v^3 ~ —2.1547 (the point where the left boundary in Fig. 5(a) 
crosses the real u axis) and apparent divergences at u = 1 and u = — 1; however, none 
of these singularities actually occurs since all of these points are outside of the complex- 
temperature extension of the FM phase where the expression ([4.1 5|) holds. This is clear from 
the phase diagram in Fig. 5(a). 

The characteristics of the points where M vanishes continuously are listed in Table 2, 
which includes a comparison with the three 2D homopolygonal lattices. The column marked 
z gives the symbol, if one was assigned, for each of the zeros (or divergences) in the uniform 
(and, where relevant, staggered) magnetisation. The column marked "value" lists either 
the analytic expression or, where this is too long to fit into the space, a reference to the 
equation(s) where this is given in the text. The column marked adj. phases lists the phases 
which are adjacent to the given point. For conciseness, some lines list a point and its complex 
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conjugate together; in these cases, the notation 0^*^ in the adjacent phase column means 
O for the first point and O* for its complex conjugate. One may recall that the divergence 
in M for the triangular lattice occurs at z = ±l/v3 which are the endpoints of two line 
segments protruding into the complex-temperature FM phase; this is indicated in the table 
by the notation "protrusion in FM". Other notation is explained in the table caption. It is 
also of interest to perform this comparison for the lattices with even coordination number, 
where all quantities can be expressed solely in terms of u = z 2 ; this is done in Table 3. 



5 (3 • 12 2 ) Lattice 

The (zero-field) free energy for the (3-12 2 ) lattice with isotropic couplings was given implicitly 
in [^3], [I7| . For various generalisations to unequal couplings, the free energy was given in 
Refs. [fM ] and pS] . Evaluating some intermediate expressions and performing some algebra, 
we obtain the explicit expression, in terms of z (and K = —(1/2) Inz), 

h-12 = \*C+\ ln (! + *) 



i r r d6ido 2 , 

+ tt: / "t^tt In 



^3-12 + -B3-12-P( 



(5.1) 



12 J-*J-« (2tt) 5 

where 

A 3 _ 12 = 1-Q Z + 24z 2 - 42z 3 + 66/ - 42z 5 + 48/ - 6z 7 + 21z s (5.2) 

B3-12 = -2z(l -z + 2z 2 )(l + z)(l - zf (5.3) 

Since this lattice has odd coordination number q = 3, Theorem 6 states that, in addition 
to the trivial singularities at K = ±00, i.e., z = 0, 00, the free energy / has a singularity at 
z = —1. This is evident in eq. ( |5.1| ). If one takes the branch cut for the logarithm to extend 
from z = — 1 to z = — 00, then as one approaches z — — 1 from the direction of the origin, / 
becomes negatively infinite. In the case of the honeycomb lattice, as was noted in Ref. ||13|| , 
this point lies on the continuous locus of points where the argument of the logarithm in the 
integrand analogous to ( p.l| ) vanishes. However, in contrast to the honeycomb lattice, for 
the (3 • 12 2 ) and (4 • 8 2 ) lattices, this singularity is an isolated one. Aside from these isolated 
singular points, the free energy also has non-analyticities which arise from the vanishing of 
the argument of the logarithm inside of the integral in eq. ( |5.1| ). The condition for this 
vanishing is 

^3-12 + 5 3 -i 2 x = ( 5 - 4 ) 
where x represents P(9\, 62) as given above in eq. (|4.2|) . This yields the curve shown in Fig. 
6(a). As marked there, these curves serve to bound the complex-temperature extensions of 
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Table 2: Points at which M and M st vanish continuously (or diverge) for the (3 ■ 6 ■ 3 • 6) 
(kagome), (3-12 2 ), and (4-8 2 ) (bathroom tile) lattices, together with the three homopolygonal 
2D lattices, for comparison. Under the M and M st columns, an entry marked — means that 
the point cannot be reached from within the FM and AFM phases, respectively. The symbol 
N means that the model has no AFM phase, (mp) means a multiple (intersection) point 
through two or more phase boundary arcs pajs^. See text for further discussion. 
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Table 3: Points in the u plane at which M and M st vanish continuously (or diverge) for the 
lattices with even coordination number q. Notation is the same as in Table 2. 

the FM and PM phases and also two phases which have no overlap with any physical phase, 

and are labelled O and O*. The border of the complex FM phase crosses the real z axis at 
the points 

z a ± = \ [-(V3 + 1) ± (4 + H V3)V2] (5.5) 

with numerical values 

z a+ = 0.19710468... (5.6) 

z a _ = -2.9291555... (5.7) 

The physical critical point of the model is z c ,3-i2 = z a+ . The intersection points nearest to 
the real z axis are given by z r and z*, where 

1 - 2i 



5-^6^ (5i 



where 8^ was given above in eq. ( |4.7|) (for certain formally analogous points in the u plane). 
The insection points farther out from the real z axis occur at z = ±z. For x = 3, there are 
four double roots at z a+ , z a _, and z|, where 



1 

Zb =2 



(v / 3-l) + (4-fv / 3) 1/2 " 

0.36602540 + 0.66586461i (5.9) 
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As can be seen from Fig. 6(a), the points z& and zl lie roughly at the centers of the 
outer boundaries of the respective O and O* phases. As x decreases from 3 in the range 
— 1 < x < 3, these double roots split apart into pairs the members of which move away 
from these four points, tracing the curves shown in Fig. 6(a). At x = — 1, these points 
rejoin at the four intersection points z = ±z and z = z r and z*. As x decreases further 
in the range —3/2 < xj < —1, these four double roots split again, with two of the roots 
moving toward each other on the arc forming the inner boundary of the O phase adjacent 
to the FM phase, and similarly with two other roots on the inner boundary of the O* phase, 
while the other four roots move outward away from the origin along the two arcs and their 
complex conjugates. Finally, at x — —3/2, two pairs of roots on the inner arcs join to form 
double roots at z = ±i/\/3, while the four outer roots reach the endpoints of the arcs. These 
endpoints are as follows: 



z e± = - 



1 W3 ± ( e 2W3 _ (4/3)^-/6 



1/2' 



(5.10) 



with the values 



z e+ ~ 0.29556791 - 0.3588689i (5.11) 

z e _ ~ 0.20443209 + 1.2248943? (5.12) 

are, respectively, the endpoints of the upper parts of the arcs in the upper and lower half z 
plane. Their complex conjugates z* + and z*_ form the endpoints of the lower parts of these 
arcs. 

From implicit expressions given in Ref. [23|, we can express the spontaneous magnetisa- 
tion as 

(1 - 6z + 6z 2 - 6z 3 - 3z 4 )^ s (l + 3z 2 ) 1/4 (l + z 2 ) 3/s (l - 2z + hz 2 )^ 
M3 " 12 " (l-*)(l + *)V2(i_ 3 + 2*2) (5 ' 13) 

By analytic continuation, this expression holds throughout the complex-temperature ex- 
tension of the FM phase, which we have determined, as shown in Fig. 6(a). Note the 
factorization 

1 - 6z + 6z 2 - 6z 3 - 3z 4 = -3(z - z a+ )(z - z a -){z - z b )(z - z* b ) (5.14) 

Moreover, similar to another formula above (for u rather than z), we have the factorization 
1 - 2z + 5z 2 = (1 - z/z r )(l - z/z*). Finally, we note that 1 - z + 2z 2 = 2(z - z )(z - z*), 
where 



z, 



o 



-(1 + iV7) = 0.25 + i0.66143783... (5.15) 



4 



19 



These points z Q and z* lie in the middle of the two complex conjugate O phases. Thus, 
M 3 „ 12 vanishes continuously at the eight points z a+ , z a _, ±i/V3, ±i, z r , and z*. In Table 
2 we list these zeros, together with the corresponding exponents (3 and the phases which 
are adjacent at each point. Elsewhere along the boundary of the complex-temperature FM 
phase, M vanishes discontinuously. M^\ 2 also has a divergence at at the isolated point 
z = — 1 in the interior of the complex-temperature FM phase, where the free energy itself 
is also singular. In addition to these actual singularities, the formal expression (|5.13 ) has 



apparent zeros at the points z b _ = z£ + , z a , and z*, but these are not zeros of the actual 
magnetisation, since these points cannot be reached from within the complex-temperature 
FM phase where the expression ( |5.13| ) holds. As one can see from Fig. 6(a), z& lies on the 
boundary separating the O and PM phase, and similarly z% lies on the boundary between 
the O* and PM phase. The points z Q and z* lie in the interior of the O and O* phases, 
respectively. 

6 (4 • 8 2 ) Lattice 

In contrast to the (3 ■ 6 • 3 ■ 6) and (3 • 12 2 ) lattices, the (4 • 8 2 ) (bathroom tile) lattice 
is bipartite, and the physical phase diagram consists of a PM, FM, and also AFM phase. 
One way to obtain the (zero-field) free energy is from the free energy for the dual [4 • 8 2 ] 
(union jack) lattice. An explicit calculation of the latter for the case of equal spin-spin 
couplings was given in p6[ . (See also the remarks in Ref. |16[]; for the union jack lattice 



with general couplings, which we shall not need here, see also Refs. [£7]] and p8| .) Using 
this duality relation, one easily obtains the free energy for the (4 • 8 2 ) lattice, in terms of z 
(and K — —(1/2) In z), as follows: 

/ 4 _ 8 = \K + \ ln(l + z) 

+ 0/ / /o \? ln A A _ 8 (z) + B 1 . A _ 8 (z)(cos6 1 + cos6 2 ) + B 2 . 4 _ 8 (z)cos8 1 cos8 2 (6.1) 
8 J-n J-7T {2n) z L J 

where 

A 4 _ 8 (z) = (1 + z 2 ) 2 (l -Az+ 10z 2 - 4z 3 + z 4 ) (6.2) 

£i;4- 8 (*) = 2z(l - zf{\ + z){\ + z 2 ) (6.3) 

B 2 , 4 ^{z) = -Az 2 {l-zY (6.4) 

In addition to the trivial singularities at K = ±oo, i.e., z = 0, oo, and the additional isolated 
singularity at z = —1, the free energy has non-analyticities where the argument of the 
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logarithm in the integrand of eq. (|6.1| ) vanishes. These are given by the equation 



Astsiz) + Bi-4_ 8 (z)(cos9i + COS 62) + B2-4-s(z) COS^i COS 62 = 



(6.5) 



We find that the locus of points which are solutions to this equation fall, in general, into 
areas rather than curves, in the z or, equivalently, the v planes. These areas degenerate 
into points at certain special locations. In Fig. 7(a) we show the solution to eq. (|6.5|) in 
the z plane, for a grid of values of 9\ and 82- Because of the finite grid of values of Q\ 
and 62 used to make the plots, the zeros show a striped structure in certain regions; it is 
understood that in the limit where the aforementioned grid of values of 9\ and 82 becomes 
infinitely fine, these stripes would merge into coherent areas; the boundaries of these areas 
are easily inferred from the plot. The phases include the complex-temperature extensions of 
the physical FM, PM, and AFM phases, as marked, together with five phases which have no 
overlap with any physical phases, labelled Oi, O2, O2, O3, and O3. The rest of the z plane 
involves various continuous regions of points where the free energy is non-analytic (forming 
the thermodynamic limits of zeros of the partition functions for finite lattices). We shall 
discuss the special locations where the areas of zeros of the partition function degenerate 
into points below, in conjunction with an analysis of the spontaneous uniform and staggered 
magnetisations. The corresponding diagram in the v plane is shown in Fig. 7(b). As in Fig. 
7(a), because of the finite grid of values of 9\ and 62 used to make the plot, the zeros exhibit 
a striped or dotted structure in certain regions; again, it is understood that for a dense set 
of values of 8± and 62, the zeros in these regions form a continuous set. 

To our knowledge, this finding constitutes the first known case of an Ising model with 
isotropic spin-spin couplings, where the non-analyticities of the free energy (aside from the 
trivial ones at K = ±00 and, for odd q, at z = — 1, as proved in Theorem 6) form a two- 
dimensional, rather than one- dimensional, algebraic variety, i.e., lie in areas rather than on 
curves (including line segments). It is easy to understand the basis for this result. The locus 
of points at which the free energy is non-analytic (aside from the above-mentioned trivial 
singularities) is determined by the condition that in the integral occuring in the free energy, 
the expression in the argument of the logarithm vanishes. For homopolygonal lattices, this 
expression reduces to the equation 



for isotropic spin-spin couplings, where x represents the function P(#i,#2) = cos#i + cos 82 
for the square lattice and P(#i,#2) = cos^i + cos 62 + cos(6 l i + 9 2 ) for the triangular and 
hexagonal (honeycomb) lattices. (Thus for the square lattice, —2 < x < 2, while for the 




(6.6) 
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triangular and hexagonal lattices, —3/2 < x < 3.) The solutions to eq. ( |6.6|) form a one- 
dimensional algebraic variety, specifically, algebraic curves (including possible line segments). 
In contrast, for unequal spin-spin exchange constants Jj, where i = 1, 2 for the square lattice, 
and i = 1,2,3 for the triangular and hexagonal lattices, the condition for the argument of 
the logarithm in the integrand to vanish is of the form 

A sq + B sq] i cos 6x + B sq . 2 cos 6 2 = (6.7) 

for the square lattice, and 

A A + B A]1 cos #i + B A . 2 cos 6 2 + B A]3 cos(#i + 2 ) = (6.8) 

for the triangular and hexagonal lattices, where in each case, the various A and B functions 
depend on the z iy with Zi = e~ 2Ki and Ki = [3Ji. Now let the (anisotropic) ratios of the Jj's 
be fixed. Then eqs. ( |6.7|) and ( |6.8|) depend on two independent real (periodic) variables, 6*i 
and 6 2 ■ It follows that the solutions, in general, form a two-dimensional algebraic variety, 
specifically areas (which may degenerate to points at special values of z). 

The heteropolygonal lattices (3 • 6 • 3 • 6) and (3 • 12 2 ) are similar to the homopolygonal 
lattices in this respect, i.e., if the spin-spin couplings on each link are equal, then the condition 
for the argument of the logarithm in the integrand to vanish is of the form (|6.6|) . However, 
as is evident from eq. (|6.1|) , this is not the case for the (4 ■ 8 2 ) lattice. That is, even if the 
spin-spin couplings are equal for each link, the condition for the vanishing of the argument 
of the logarithm in eq. ( |6.1|) is of the form (|6.5| ), which depends on two independent real 
(periodic) variables 6\ and 9 2 . 



The expression for the spontaneous magnetisation was conjectured by Lin et al. [22] and 



proved by Baxter and Choy | 27fl . As with the other heteropolygonal lattices A, M can be 



written (where it is nonzero) as a prefactor times (1 — (/c <5 a) 2 ) 1//8 - Let us define (using the 
subscript 4 - 8 to denote (4 ■ 8 2 )) 

8z 2 (l-2z + 4z 2 -2z 3 + z 4 ) 

k <M-8) = 7j , ( 6 - 9 ) 



and 



1-z 4 (l + z 



[l + z 2 ) 



Then M 4 _ 8 = a(l - (A;<) 2 ) , i.e., 



{l + z) 1 / 2 {l-z + 3z 2 + z' i ) 1 / 2 

1/8 



(6.10) 



(1 - Az - z 4 ) 1/8 (l + 4^ 3 - z 4 ) 1/8 (l + ^ 2 ) 1/2 (1 - 2z + 6z 2 - 2z 3 + z A ) 1 ' 4 

4-8 ~ (l-z)(l + z) 1 /2(l- Z + 3z 2 + z3)l/2 t 6 " > 
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within the physical FM phase and, by analytic continuation, throughout the complex- 
temperature extension of the FM phase; M 4 _ 8 vanishes identically elsewhere. We observe 
first that under the transformation z — > 1/z, the numerator of this expression goes into itself 
times a power of z. It follows that the set of formal zeros of this numerator is invariant 
under the mapping z — > 1/z. In particular, under this mapping, the first and second quartic 
expressions in the numerator of M4_g are interchanged, up to an overall power of z: z 1/z 
=>- (1— 4z — z 4 ) — > —z^ 4 (l+Az 3 — z 4 ). The roots of these two quartics are therefore reciprocals 
of each other. Under this mapping, the other two factors, (1 + z 2 ) and (l — 2z + 6z 2 —2z 3 + z 4 ) 
transform into themselves, up to overall powers of z, and consequently, their zeros also come 
in reciprocal pairs (since, furthermore they do not have either of the self-reciprocal numbers 
z = ±1 as roots). We note the factorizations 



Az-z 4 



1 - V2(V2 + l)z - (V2 + l)z 2 } \l - V2(V2 - l)z + (V2 - 1 



- :! z 4 



1 + Az 
1 - 2z + 6z 2 - 2z 3 + z 4 
The expression ( |6.12|) has zeros at 

zi± = - 



1 + V2z -(V2- l)z 2 } \l - V2z + (V2 + 1) 



I 1 " 



2e iw/3 z + z< 



2e- in/3 z + z 2 



Z2± ~2 



while the expression ( |6 . 13|) has zeros at 



-i 



.-1 



V2± (-2 + 4v / 2) 1/2 
v / 2±i(2 + 4 v / 2) 1/2 

2 + V2± (lO + Sv 7 ^) 1 



2- v / 2±(10-8\/2) 



/2 



1/2 



The expression ( |6.14 ) has zeros at z^± and z%±, where 



e in/3 ± _ x y 



3 2ivr/3 _ ]\l/2 

Note that z$- = z% + . These roots have the approximate numerical values 

z 1+ = 0.2490384 , = 4.0154454 

zi_ = -1.663252 , Ozi-)- 1 = -0.60123183 
z 2± = 0.7071068 ± 1.383551i 



23 




-1 



0.2928932 ± 0.5730856i 



(6.23) 



^3+ 



0.8406250 + 2. 137255i 



(6.24) 



= [z 3+ 



0.1593750 - 0.4052045i 



(6.25) 



Only a subset of these formal zeros are zeros of the actual magnetisation, namely those which 
can be reached from within the complex-temperature FM phase (they lie on the border of 



points which cannot be reached in this manner, where consequently the formula does not 
apply. The true zeros of M are at z c = z± + , l/^i-, and ;?3_. The first two of these 
are, respectively, the physical critical point and the point at which the left boundary of the 
complex-temperature FM phase crosses the real z axis. As one can see from Fig. 7(a), z%_ lies 
where certain areas of zeros of the partition function degenerate to a single point separating 
the complex-temperature FM and O2 phases, and similarly for the complex conjugate point 
2:3- separating the FM and phases. Indeed, these four points are precisely the full set 
of points on the border of the complex-temperature FM phase where various areas of zeros 
degenerate to single points. The exponents with which M4_§ vanishes continuously at these 
zeros are all equal to 1/8. As one moves upward, away from the origin along the Im(z) 
(vertical) axis in Fig. 7(a), one encounters a boundary followed by a dense set of points 
where the free energy is non-analytic before one reaches the point z = i. This boundary 
prevents analytic continuation, so that the zero at z = i, and similarly, the zero at z = — i, 
in the expression ( |6.1 1| ) are not true zeros of M. Note that the points ±i lie where the Oi 
and PM phases are directly adjacent, and separated by these points. Similarly, all of the 
remaining zeros in the numerator of the expression ( |6.11| ) are not actual zeros of M. However, 
it is of some interest to observe where these spurious zeros lie. The point z<i+ is located where 
the border between the O3 and PM phases narrows to zero thickness (elsewhere this border 
is comprised of arcs forming strips, as one can see from Fig. 7(a)). An analogous comment 
applies, mutatis mutandis, for the complex conjugate point z% + and the border between the 
O3 and PM phase. The reciprocal point I/22- is located where the border between the O2 
and PM phase narrows to zero thickness, and similarly for the point l/z% + and the border 
between the and PM phases. The point z% + lies where two areas of zeros of the partition 
function degenerate to a point separating the O3 and AFM phases, and similarly for z\ and 
the border between the O3 and AFM phases. Finally, none of the the apparent divergences 
in M 4 _ 8 actually occurs; z = 1 lies in the interior of the PM phase; z = — 1 is in the interior 
of the Oi phase, and among the roots of the cubic 1 — z + 3z 2 + z 3 , z ~ —3.38 is in the AFM 
phase, while z ~ 0.191 ± 0.509z lie in the interiors of the O2 and O2 phases, respectively. 
In the case of the point z = — 1, it is interesting to compare this with the situation for the 



this phase) where the formula ( p. 11 ) holds; the other zeros are spurious, since they occur at 
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(3 • 12 2 ) lattice. In both cases, since q is odd, / is singular at z = — 1. For the (3 • 12 2 ) lattice, 
this point lies within the complex-temperature extension of the FM phase, and M is also 
singular (divergent) at this point. In contrast, for the (4 • 8 2 ) lattice, this point is not within 
the complex-temperature FM phase, but rather in the interior of one of the O phases, and 
hence there is no associated singularity in M at this point. 

As with other bipartite lattices, the staggered spontaneous magnetisation M sti4 _ 8 is given 
by eq. in terms of M^-g in the complex-temperature extension of the AFM phase and 

vanishes elsewhere. We may thus read off the zeros of M s4 4 _ 8 immediately; these are precisely 
the inverses of the four zeros of M, i.e., z~ l , (z^)^ 1 = z 3+ , and (z^) -1 = z$ + . The first 
point is the physical critical point separating the PM and AFM phase. The second is the 
point where the boundary between the AFM and O phases crosses the negative real z axis. 
The third is where the border between the O3 and AFM phases narrows to zero thickness, 
and similarly for the fourth point, separating the O3 and AFM phases. The actual zeros of 
M and M st are listed in Table 2. 

In Table 4 we compare the general properties of the complex-temperature phase diagrams 
for the Ising model on homopolygonal and heteropolygonal 2D Archimedean lattices, in 
terms of the variable z, or equivalently, v. Since we have already indicated in Table 2 which 
lattices have AFM phases, and since all have PM and FM phases, we do not include that 
information here. The entries in the column marked No are the number of O phases. We 
observe that the Ising model on heteropolygonal lattices exhibits more of these O phases 



than on homopolygonal lattices. As we have discussed in our earlier papers [11, O, O 



the points at which the curves (or line segments, if present) cross are singular points of the 
algebraic curves, in the technical terminology of algebraic geometry [21|. This usage of the 
term should, of course, not be confused with the different meaning of "singular point" in 
statistical mechanics; for example, the point z c = (—1 + 2/v / 3) 1 ^ 2 of the Ising model on 
the kagome lattice is a singular point (actually, the physical critical point separating the 
PM and FM phases) in the statistical mechanical sense, but is not a singular point of the 
curve forming the outer boundary of the complex-temperature FM phase, in the algebraic 
geometry sense (see Fig. 5(b)). A singular point of the algebraic curve is denoted as a 
multiple point of index n& if rib branches (arcs) of the curve pass through the point. (We 
use the terms "multiple point" and "intersection point" synonymously here.) In the column 
marked "mult, pt.", we list information about these multiple points, including the number 
and the index. The symbol cc means that the multiple points occur in complex conjugate 
pairs, while "real" means that the multiple points occur on the real axis in the z (equivalently, 
v) plane. For a normal double point, i.e., a multiple point of index 2, one may assign an 
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angle 6 cr as the angle between the tangents of the two branches (arcs) of the curve which 
pass through the point. All three homopolygonal lattices have multiple points on the curves 
(including line segments) across which the free energy is non-analytic. Each of these multiple 
points has index 2 and 8 cr = tt/2, i.e., the two branches cross in an orthogonal manner. Our 
results as shown in Figs. 5 and 6 show that the multiple points for the (3 • 6 • 3 • 6) and 
(3- 12 2 ) lattices also have index 2 and 8 cr = tt/2. Note that because the transformation (|2.3| ) 
relating z and v and also the transformation u = z 2 are conformal mappings, they preserve 
angles, so that a given 8 cr is the same for a multiple point in the z, v, and u plane. For all 
of these four lattices, the multiple points occur as complex conjugate pairs as functions of 
z. A general observation is that the heteropolygonal lattices have more multiple points on 
their curves than the homopolygonal lattices. The situation with the (4 • 8 2 ) lattice is the 
most complicated, since in this case the loci of points where the free energy is non-analytic 
form areas instead of curves. As one can see from Fig. 7, there are a number of multiple 
points through which more than one of the boundaries of these areas pass. There are 18 
such points in all; of these, 14 consist of complex conjugate pairs while the remaining four lie 
on the real z or v axis. Of the 14 complex conjugate multiple points, 12 have index rib = 2 
(with various values of crossing angles 6 cr ). The four multiple points lying on the real axis 
(namely, z c = zi+, l/z c , zi_, and l/z\-) have index rif, = 2 and 8 cr = 0. In the terminology of 
algebraic geometry, they are therefore (simple) tacnodes. We recall that a (simple) tacnode 
is defined as a multiple point of an algebraic curve through which two branches pass, in an 
osculating manner, i.e. such that their tangents coincide where the branches touch, so that 
rib = 2 but the number of distinct tangents at the point is n t = 1 [^lj . We have denoted 
this in Tables 2 and 4 by the notation "tn" , standing for "tacnode" . A generalised tacnode 
is a multiple point of an algebraic curve for which some of the tangents of the branches 
passing through the point coincide, so that rib > n t (and n b may be greater than 2). The 
(4 • 8 2 ) lattice is the first one which we have studied which exhibits tacnodal points on the 
loci of points where the free energy is non-analytic. We come next to the multiple points at 
±i. These are, again, of a type unprecedented in any of the cases which we have previously 
studied; they have index 4 and are the first example of a multiple point of index higher than 
2. Of the 4 branches which pass through this point, two cross at an angle of tt/2 with respect 
to each other, the northeast - southwest and northwest - southeast curves), while two others 
pass through the points in a north - south direction, hence with a crossing angle of it /A with 
respect to the former two curves and a crossing angle of with respect to each other. This 
multiple point is therefore a higher-order tacnodal point, with = 4 and n t = 3. Finally, 
in the last column of Table 4 we list the number of endpoints of the algebraic curves and in 
which complex-temperature phase these endpoints lie. For all cases, the endpoints occur in 
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lattice 


AT 


mult. pt. 


endpt . 


sq (4 4 ) 


1 


2: 2, cc 





tn (3 b ) 


1 


2: 2, cc 


2, cc, in FM 


nex ) 


n 
U 


Z. z, cc 


z, CC, 111 JrlVl 


(3 • 6 • 3 • 6) 


2 


4: 2, cc 


4, cc, in PM 


(3 • 12 2 ) 


2 


4: 2, cc 


4, cc, in PM 


(4 • 8 2 ) 


3 


18; see below: 









12: 2, cc 








2: 3, cc, tn (±i) 








4: 2, real, tn 





Table 4: Comparison of some properties of complex-temperature phase diagrams, as func- 
tions of z or equivalently v, for the Ising model on homopolygonal and heteropolygonal 
Archimedean lattices. No denotes the number of O phases. In column marked "mult, pt.", 
the entry 2: 2, cc means that there are 2 multiple points, each of index 2, and they occur as 
complex conjugate pairs, etc. for other entries. The notation "tn" means a tacnodal multiple 
point. For the (4 • 8 2 ) lattice we have listed the properties of the various types of multiple 
points underneath the general number, 18. See text for further discussion. 

complex conjugate pairs. 

For lattices with even q, viz., the square, hexagonal, and kagome lattices, a more compact 
way to present the complex-temperature phase diagrams is in the u plane. In Table 5 we list 
the analogous characteristics for these lattices. 

7 Conclusions 

In this paper, using exact results, we have determined the complex-temperature phase dia- 
grams and singularities of the spontaneous magnetisation on three heteropolygonal Archimedean 
lattices, (3 ■ 6 ■ 3 • 6) (kagome), (3 • 12 2 ), and (4 • 8 2 ) (bathroom tile). This study reveals a rich 
variety of complex-temperature singularities and provides an interesting comparison with 
the situation for the three homopolygonal lattices. In particular, we have found the first 
example of a lattice where, even for equal spin-spin exchange couplings, the nontrivial non- 
analyticities of the free energy lie in areas rather than on curves, in the z (or equivalently, 
v) plane. We have also given a simple explanation of why this happens. 
This research was supported in part by the NSF grant PHY-93-09888. 
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lattice 


No 


mult. pt. 


endpt. 


sq (4 4 ) 





1: 2, real 





tri (3 6 ) 





1: 2, real 


1, real, in FM 


(3 ■ 6 ■ 3 • 6) 


1 


2: 2, cc 


2, cc, in PM 



Table 5: Comparison of some properties of complex-temperature phase diagrams, as func- 
tions of u, for the Ising model on lattices with even coordination number q. Notation is as 
in Table 4. 
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Figure Captions 

Fig. 1. The (3 • 6 • 3 • 6) (kagome) lattice. 

Fig. 2. The (3 • 12 2 ) lattice. 

Fig. 3. The (4 • 8 2 ) (bathroom tile lattice. 

Fig. 4. The [4-8 2 ] (union jack) lattice, which is the Laves lattice dual to the Archimedean 
lattice (4-8 2 ). 

Fig. 5. Complex-temperature phase diagram for the Ising model on the 3 ■ 6 ■ 3 ■ 6 
(kagome) lattice, in the variable (a) u (b) z (c) v. 

Fig. 6. Complex-temperature phase diagram for the Ising model on the 3 • 12 2 lattice, 
in the variable (a) z (b) v. 

Fig. 7. Complex-temperature phase diagram for the Ising model on the 4 • 8 2 (bathroom 
tile) lattice, in the variable (a) z (b) v. 
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